Formation of hyperfine fields in alloys 
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This work deals with the analysis of experimental data on the average magnetization of Fe\- x Me x 
(Me=Sn,Si) disordered alloys, the average and local hyperfine fields (HFF) at the Fe nuclei. The 
effect of the metalloid concentration on the HFF is studied with the help of the results of first- 
principles calculations of ordered alloys. The disorder is taken into account by means of model 
systems. The dependences obtained correspond to those experimentally observed. Experimental 
data on the ratio of the average HFF at Fe nuclei to the average magnetisation in alloys with 
sp-elements show that the ratio decreases proportionally with the metalloid concentration. This 
change in the ratio is bound up with three factors. First, the contribution of the valence electron 
polarization by the neighboring atoms, that is positive (unlike the polarization by the own magnetic 
moment), increases with the change of the disorder degree (increase of concentration). Second, the 
appearance of the impurities, i.e. metalloid atoms, in the nearest environment of Fe leads to the 
orbital moment increase. And, finally, the change of the disorder degree, as in the first case, results 
in an increase in the orbital magnetic moment and its positive contribution to the HFF. The value 
and the degree of the influence of these contributions to the HFF is discussed. 
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In the experimental physics of solids there are few 
methods for measuring local characteristics. Prominent 
among them are nuclear techniques (e.g. Mossbauer 
spectroscopy) determining the hyperfine magnetic fields 
(HFF) at nuclei. It is widely believed that these tech- 
niques can reflect some chemical and topological features 
of the atomic surroundings of the excited nuclei. In- 
deed, in numerous studies these spectra are interpreted in 
terms of phenomenological models considering only the 
nearest atomic environment (see e.g. [1-3]). There is 
no doubt that such a description is often useful and effi- 
cient. We think, however, that there are some cases when 
such a simplified approach is not justified. Moreover, for 
the spectra of transition metals and alloys on their ba- 
sis, the successful description within the framework of 
such models is rather an exception than a rule, and ev- 
ery time additional arguments are needed to justify these 
limitations. For example, such an approach makes some 
sense in the case of disordered systems due to strong lo- 
calization of the d- electrons responsible for the itinerant 
magnetism of transition metals [4]. This is one of the 
reasons why sometimes the spectra of disordered alloys 
are successfully described by the Jaccarino- Walker- type 
models where the local magnetic moments (LMM) and 
HFFs are assumed to be proportional to the number of 
metalloid atoms in the nearest environment. 

In our opinion, even in this case the experimental spec- 
tra may provide more reliable and complete information 
without the use of phenomenological models and restric- 
tion to the nearest environment. One should then analyze 
the HFF peculiarities from the " first- principles" calcula- 
tions. This work is devoted to such an analysis of the 
hyperfine magnetic fields, magnetization and local mag- 
netic moments and their interrelation in the disordered 
Fe\- x Sn x , Fe±- x Si x alloys. Here the magnetic charac- 
teristics of clusters with a given impurity configuration in 



disordered alloys are taken from the first-principles calcu- 
lations of the translationally invariant systems, supposing 
the interaction between the clusters to be not significant 
because of small free length of electron in these alloys. 
Nowadays, there is no possibility to conduct these cal- 
culations for the disordered systems, moreover, it should 
be noted that even calculations of the ordered alloys of- 
ten do not give the required agreement with experiment 
and reveal only the main features of the HFF and LMM 
behavior. 

The calculations were performed by the full-potential 
linearized augmented plane wave method (FLAPW) us- 
ing the WIEN-97 program package [5]. The results are 
presented in Table 1. 

The systems were simulated on a BCC lattice which 
in the disordered alloys under consideration is retained 
within a wide concentration range [1]. The lattice pa- 
rameters were chosen in accordance with the experi- 
mental values for x=3.125 at.% and x—6.25 at.%. It 
should be mentioned that even at small concentrations 
the BCC lattice is somewhat distorted due to the re- 
pulsion/attraction by the Sn/Si atom of the surround- 
ing Fe atoms. As shown in our paper [6], the changes in 
magnetic characteristics because of this relaxation are in- 
significant. Though the results in Table 1 were obtained 
with allowance for this relaxation, we do not discuss it 
here. 

1. AVERAGE AND LOCAL MAGNETIC MO- 
MENTS. 

Fig.l presents the experimental data and the calcu- 
lated averages of the magnetic moment per one Fe atom. 
The average magnetic moment M — Y2iMd l + M mt , 
where Md % is the spin magnetic moment over the muffin- 
tin (MT) sphere of the i-th Fe atom (hereinafter we will 
refer to this value as local magnetic moment, LMM), 
M mt is the spin magnetic moment over the unit cell with- 
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FIG. 1. The experimental data and the calculated averages 
of the magnetic moment per one Fe atom 

out the MT- spheres. The main contribution to the LMM 
comes from the d-electrons that are almost entirely inside 
the MT-sphere, whereas M int is formed by the s- and p- 
electrons and has a small negative value as compared to 
the LMM. We do not take into account the contribution 
of the orbital magnetic moment that is about 0.045 in 
our calculations. In experiments this value is 1.5?2 times 
larger and comprises 0.08 [is [7]. It is believed that 
the inclusion of the orbital polarization in the exchange- 
correlation potential [8] improves the agreement between 
the calculations and the experimental values, but here we 
didn't use such a potential supposing that the calculated 
value can be multiplied by the factor close to two. The 
main reason for the neglect of the orbital contribution 
lies in the fact that its variations with concentration and 
configuration of metalloid atoms are small as compared 
to M even with allowance of the actual value and amounts 
less than 1 % (see the last but one column in Table 1), 
whereas the Mdi variations range up to 15 %. Fig.l dis- 
plays a rather good agreement of the magnetic moment 
with the experimental data, though the theoretical values 
of the magnetic moment are somewhat higher than the 
experimental ones. We believe that the disorder, which is 
not taken into account here, reduces the magnetization 
by 2-3 % [9]. Let us mention some peculiarities of the 
magnetic moments formation. 

A. The magnetic moment does not depend on the met- 
alloid type and, as shown in Ref.10, is governed by the 
lattice parameter, a. At equal concentrations, the lattice 
parameter of the Sn alloy is greater than that of the Si 
alloy, so the magnetic moment in the first case is greater. 

B. The LMM of the Fe atom closest to the impurity is 
of the smallest among others value (Table 1). As noted 
in Ref.10, this difference is defined by the competition 
between two mechanisms: the LMM reduction due to 
flattening of the d-band because of the s-d hybridization 
that is the strongest near the impurity, and the LMM 
increase due to narrowing of the d-band because of a de- 
crease of the wave function overlap. There also exists 



third mechanism of the LMM reduction due to the dif- 
ference in the impurity-potential screening by d-electrons 
(the difference in pushing out the impurity levels by the 
bands with spin up and down), which was revealed in this 
work by a comparison of the number of d-electrons within 
the MT-sphere for different Fe-atom positions (there are 
more d-electrons near the impurity, Table 1). 

C. The LMMs are concentration dependent. So the 
LMM of the Fe atom closest to the Si atom is 2.181 
/i b at x=3.125 at.% (Fe 31 Si), and 2.262 /jb at x=6.25 
at.% (Fei^Si). This LMM increase with concentration 
holds in general for all the non-equivalent positions (Ta- 
ble 1). However, in spite of the LMM increase, the aver- 
age magnetic moment in the Fe\- x Si x system somewhat 
decreases 2.238 [i B -I 2.230 [i B (Fig.l). This is the re- 
sult of a higher probability of finding a Fe atom with a 
neighboring impurity, which increases the number of Fe 
atoms having lower LMM values. 

2. HYPERFINE MAGNETIC FIELDS AT NUCLEI. 

The program package WIEN-97 allows one to calculate 
the interaction between the nucleus magnetic moment 
and the spin and orbital magnetic moments of the elec- 
tron subsystem. The spin dipolc contributions are small 
( 2 v 3 kGs), they are suppressed due to the symmetry 
relations, and the main contribution comes from the spin 
polarization at the nucleus (Fermi-contact interaction) 
and the orbital magnetic moment. The Fermi-contact 
interaction may be divided into Hf ore (the core-electron 
polarization) and H^ al (the valence-electron polariza- 
tion) , and hence, the resulted field at the i-th site contains 
three terms: H l = Hf ore + H™ 1 + H° rb . For the core- 
electron polarization the simple relation Hf ore = j s Mdi 
is satisfied, where -f s does not depend to a high accuracy 
on neither the metalloid type nor its concentration, and is 
determined only by the approximation for the exchange- 
correlation potential. Here we use the GGA approxima- 
tion [11] which gives j s w 123fcGs/ hb- In our opinion 
just this simple expression between the LMM and Hcore 
makes it possible to use phenomenological models ne- 
glecting the effect of the atoms in the second, third, etc. 
coordination spheres. Really, from Table 1 one can see 
that the main distinctions in LMM at a certain concen- 
tration arc connected with the presence or absence of 
the metalloid atom in nearest environment. So, this also 
determines the variations in Hf ore . Though the experi- 
ments showed that the proportionality between the mag- 
netic moment and HFF is not as good as we found, and 
the proportionality coefficient essentially decreases with 
concentration. As we shall see later this is connected 
with the other two contributions to the HFF. 

The proportionality Hf rb = j° rb M? rb is fulfilled in a 
somewhat worse way, but still rather satisfactorily. Note, 
however, that j orb is positive and about five times larger 
in magnitude than 7 s . Hence, if the changes of M° rb af- 
fect the magnetic moment only slightly, the changes in 
Horbi may amount to 20 kGs even at low concentrations 
(see Table 1). The nature of these changes in disordered 
alloys is discussed in more detail in Ref.12. Here we men- 



2 



4 















Fe15Sn 














Fe31Sn 










vertical bars show distance in Fe15Sn 


- 




\ 








and in Fe31Sn - 




i ■ 





























4 5 6 7 8 9 10 11 12 13 14 

r(a.u.) 

FIG. 2. The polarizations of the valence electrons by a 
Fe-atom magnetic moment as a function of distance 
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FIG. 3. The experimental HO and HI as a function of con- 
centration 



tion only the main features. The H° rb increases along 
with the M° rb increase with concentration. The H° rb 
takes the largest value at the atom closest to the metal- 
loid atom. The increase of H° rb with allowance for the 
actual values of M° rb , that are twice as large as the cal- 
culated one, comprises 15-^20 kGs even at low concentra- 
tions as compared to that in pure Fe. On the strength of 
the qualitative character of the last statement, we can say 
about the tendency of the variations. Finally, as shown in 
[10], the orbital contribution increases also with disorder 
(that is, with concentration). 

The HV al behaves in a more complicated way. This 
is primarily associated with strong delocalization of the 
s- and p-like electrons that interfere at sites with dif- 
ferent magnetic and charge properties, and therefore 
the H? behavior cannot be in fact quantitatively pre- 
dicted. However, we succeeded in revealing some quali- 
tative regularities supported by experimental evidence. 
First of all, we analyzed the valence contribution us- 
ing the simple functional dependence of the magnetic 
moment screening in the RKKY (Rudcrman - Kittcl 
- Kasuya - Yosida) theory, as it was done in Rcf.9: 
H™ 1 =A + BJ2j M^cos(2k f r/T + <j))/r 3 . Such a simpli- 
fied treatment of H™ 1 is hardly justified in our case, but 
we hope to have determined the main qualitative depen- 
dences. Solving the inverse task for the ordered cluster 
of size 200 a.u. we receive the most probable values of A, 
B, T and 4> in the alloys Fe^Sn and Fe 3 iSn. The cor- 
responding functions of Bcos(2kjr/T + (f>)/r 3 are shown 
in Fig. 2. Of special interest is the fact that for both 
concentrations the spin polarization of electrons is posi- 
tive for the I and II coordination spheres, which entirely 
contradict the results of a similar processing of the exper- 
imental data. This is due to the fact that during process- 
ing of the experimental data the difference between HO 
(HFF at the nucleus of the Fe atom without the metalloid 
atoms in its nearest environment) and HI (with one im- 
purity atom in the nearest environment) was attributed 
to the changes in H val , whereas in an alloy there are 
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FIG. 4. The configurationally averaged negative contribu- 
tion Hval 

differences in the local magnetic moment that primarily 
affect H core and H orb . Fig. 3 presents the experimental 
H a and Hi as a function of concentration. Table 1 shows 
that the difference between these quantities can be suc- 
cessfully explained by the LMM magnitude and hence 
by the core-electron polarization. However, if everything 
were determined only by Hf ore we should expect an in- 
crease in Hq and Hi with concentration in accordance 
with the LMM increase, which is not the case, as for the 
Fei- x Sn x alloy the magnetic moment increases much 
more quickly than HO and HI, and there is no increase 
at all for Si. In reality the expected increase is com- 
pensated by the decrease in magnitude of the configura- 
tionally averaged negative contribution H val . Fig. 4 gives 
the averaged values of H val of a disordered cluster of size 
200 a.u. with a certain number of impurities in the first 
coordination sphere. The averaging was performed in 
assumption that the Fe atoms are distributed randomly 
and polarize the conductivity electrons at distance r ac- 
cording to the model function for concentrations 6.25 and 
3.125 % (see Fig. 2). The decrease of the averaged Hval 
with concentration is due to the positive values of the 
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RKKY polarization at the distance of the first and sec- 
ond coordination spheres (see Fig. 2). Two main features 
in the behavior of the Hval averaged over configurations 
can be noticed. First, the magnitude of H va i at the Fe 
atom without impurities in nearest environment is less 
by 5-10 kGs than H va i at the Fe atom with one impu- 
rity in the environment, and second, the magnitude of 
the averaged H va i decreases with concentration for both 
configurations of the environment. 

Thus, the use of the "first-principles" calculations 
makes it possible to explain the main peculiarities of 
the HFF behavior in the low-concentration disordered 
alloys Fe\- x Me x . The difference between Ho and H\ 
H '- H x « -20 kGs consists of H% al - H? al w7vl0 
kGs, H$ ore - H1 ore « - 15 + -25 kGs and Hg rb - H? rb w 
—5 -. — 10 kGs. The decrease of the proportionality coef- 
ficient between the HFF and the magnetic moment with 
concentration results from the decrease of magnitude of 
H val and increase of magnitude of H orb that are opposite 
to the H core . 

We believe that the general relations obtained here 
will be useful in processing the experimental data on the 
HFFs in disordered alloys of transition metals and non- 
magnetic impurities. 

The present work was supported by Russian Founda- 
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TABLE I. The results of the calculations: Configuration of impurities in the Fe-atom environment, [nm...] denotes the 
number of metalloid atoms in the first (n), second (m) etc. spheres. Number of such Fe atoms in the unit cell, Nf e - Number of 
d-electrons in the MT sphere, N d . Magnetic moment, M d , in the MT sphere. Contribution of the core-electrons polarization 
to the HFF, H core . Contribution of the valence-electrons polarization to the HFF, H val . Orbital magnetic moment, M orb . 
Orbital contribution to the HFF, H . 





[nm...], 


N Fe 


N d 


M d , (IB 


# core ,kGs 


H val ,kGs 


M° r \ (I B 


H orb ,kGs 


Fe 31 Sn, 


[1000] 


8 


6.006 


2.246 


-277.50 


-37.59 


.0489 


27.94 


a=21.804a.u. 


[0003] 


8 


5.951 


2.507 


-310.35 


-30.94 


.0483 


26.41 




[0100] 


6 


5.964 


2.432 


-300.17 


-34.84 


.0482 


26.62 




[0020] 


6 


5.963 


2.425 


-299.79 


-22.39 


.0479 


26.28 




[0000] 


2 


5.987 


2.391 


-295.13 


-34.34 


.0441 


23.66 




[0000] 


1 


5.959 


2.423 


-299.68 


-24.17 


.0470 


25.15 


Fei 5 Sn, 


[1003] 


8 


6.114 


2.442 


-301.01 


-20.25 


.0548 


29.77 


2a=21.985a.u. 


[0200] 


3 


6.090 


2.530 


-312.64 


-29.97 


.0535 


29.08 




[0040] 


3 


6.084 


2.542 


-314.55 


-2.46 


.0550 


29.38 




[0000] 


1 


6.142 


2.390 


-294.33 


-29.97 


.0450 


23.32 


Fe 3 iSi, 


[1000] 


8 


5.999 


2.181 


-269.44 


-38.68 


.0497 


28.27 


a=21.604 a.u. 


[0003] 


8 


5.967 


2.413 


-297.89 


-35.15 


.0471 


25.55 




[0100] 


6 


5.977 


2.343 


-288.51 


-31.31 


.0446 


24.46 




[0020] 


6 


5.965 


2.385 


-294.06 


-16.92 


.0468 


25.32 




[0000] 


2 


5.965 


2.489 


-308.13 


-19.72 


.0485 


26.44 




[0000] 


1 


5.969 


2.389 


-294.81 


-13.56 


.0451 


24.28 


Fe 15 Si, 


[1003] 


8 


6.130 


2.262 


-277.20 


-29.11 


.0594 


32.38 


2a=21.585a.u. 


[0200] 


3 


6.115 


2.353 


-289.52 


-21.64 


.0456 


24.44 




[0040] 


3 


6.094 


2.432 


-298.87 


5.49 


.0509 


27.31 




[0000] 


1 


6.094 


2.536 


-313.82 


-13.89 


.0496 


26.83 



•5 



